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We report the results of an HI blind survey of 80 deg- of the Virgo cluster, based on the 08" < S < 16° strip of ALFALFA, the 
Arecibo Legacy Fast ALFA Survey. 187 HI sources of high significance are found providing a complete census of HI sources in this 
region of the Virgo cluster (-1000 < cz < 3000 km s" 1 ) with M H i > 10 75 ~ 8 M Q . 156/187 (83%) sources are identified with optical 
S_i ■ galaxies from the Virgo Cluster Catalogue (Binggeli et al. 1985), all but 8 with late-type galaxies. Ten sources are not associated with 

' optical galaxies and were found to correspond to tidally-disrupted systems (see Kent et al. 2007 and Haynes, Giovanelli and Kent 

2007). The remaining 21 (11%) are associated with galaxies that are not listed in the Virgo Cluster Catalogue. For all sources with an 
optical counterpart in the Sloan Digital Sky Survey, we analyzed i-band SDSS plates to measure optical structural parameters. We find 
that in the Virgo cluster: i) HI inhabits galaxies that are structurally similar to ordinary late-type galaxies; ii) their HI content can be 
predicted from their optical luminosity; iii) low surface brightness galaxies have low optical luminosity and contain small quantities 
of neutral hydrogen; iv) low surface brightness, massive Malinl type galaxies are comfortably rare objects (less than 0.5 %); v) there 
are no "dark-galaxies" with HI masses M H i > 10 7 ' 5 ~ 8 M Q ; vi) less than 1% of early-type galaxies contain neutral hydrogen with 
(N . M Hl > 10 7 5 ~ 8 M (di Serego Alighieri et al. 2007). 

\o . 
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ryQ ' 1. Introduction Centaurus, carried out with the Parkes multibeam system, sen- 

. sitive to Mhi > 1O 8 M , with deep optical follow-up plates. A 

L) • The faint-end slope of the halo mass function either predicted more recent one that combines a large area with a high sensi- 

.£h '. analytically (Press & Schechter 1974) (a = -1.8) or by nu- tivity is the 86 deg 2 survey in the direction of Canes Venatici 

^ . merical CDM simulations (e.g. Jenkins et al. 2001)(or = -2) (Cvn) carried out with the WSRT (Westerbork Synthesis Radio 

$— i ' is uncomfortably steeper than the observed slope of the opti- Telescope) by Kovac (2007). Incidentally the slope of the faint- 

5^ ! cal luminosity function of galaxies (e.g. Blanton et al. 2003) end mass function derived using 70 sources with M HI > 10 6 5 M o 

(a = -1.1) (see however a = -1.5 in r band by Blanton et from this survey is flatter than ever (a = -1.17). 

al. 2005). If many sterile halos (i.e. unable to give birth to stars), The Arecibo telescope, equipped with the new 7 beam system 

the so called dark galaxies, or many low mass galaxies with low ALFA, has recently started a number of HI blind surveys. The 

optical surface brightness, yet retaining some neutral hydrogen Arecibo Galaxy Environment survey (AGES) (Auld et al. 2006) 

existed, outnumbering the normal galaxies, the faint-end slope of has covered some patches of sky, including 5 deg 2 of the Coma 

the HI mass function would be steeper than that of the optical lu- supercluster containing the cluster A1367, with a sensitivity of 

minosity function, perhaps reconciling the observations with the M m = 10 8 8 M o at the distance of Coma (300 s integration time 

theoretical predictions. However the measured faint-end slope per beam, rms=0.84 mJy/beam; Cortese et al. 2007). Highly sta- 

of the HI mass function, even in its most robust determination tistically significant results will be obtained from ALFALFA, the 

(Springob et al. 2005,) is significantly flatter (a = -1.24) than Arecibo Legacy Fast ALFA Survey, that covers the Arecibo sky 

the theoretical one (see also Zwaan, Briggs & Sprayberry, 2001). (7000 deg 2 ) in a 2-pass drift-scan mode (48 s integration, with a 

The Zwaan et al. (2005) HI mass function, based on HIPASS typical rms=2.1 mJy/beam at 10 km/s resolution; Giovanelli et 

(HI Parkes All-Sky Survey) data (Meyer et al. 2004), is flatter al. 2005). Besides a robust determination of the faint-end slope 

(a = -1.37) than predicted, although it suffers from statistical of the HI mass function, (based on several hundred galaxies 

limitations due to the paucity of objects (40) with M H1 < 1O 8 M with M HI < 10 7 5 M Q ) ALFALFA will provide more evidence 

and from distance uncertainties (Masters et al. 2004). To im- on whether low surface brightness galaxies, that are so elusive 

prove the determination of the HI mass function, several blind optically (Disney 1976; Disney & Phillipps 1987; Sabatini et 

HI surveys have been recently carried out or are under way. al. 2003), will change our understanding of galaxies in the lo- 

HIDEEP (Minchin et al. 2003) is a deep HI survey of 32 deg 2 in cal universe. This issue has two aspects that could be addressed 
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and solved with ALFALFA. One is to increase the number of 
dwarf galaxies surveyed to better determine the number density 
of low surface brightness galaxies that, in spite of their low op- 
tical visibility, are HI rich. (Galaxies of decreasing mass have 
an increasing gas fraction and a decreasing optical (/-band) sur- 
face brightness (Boselli et al. 2001; Warren et al. 2006)). The 
other is to confirm and corroborate the finding of Briggs (1997), 
Rosenberg et al (2002), Zwaan et al (2003), Minchin et al. (2004) 
that the frequency of Malinl-type galaxies, i.e. massive, gas rich 
objects of very low optical surface brightness, is insignificant 
(see also Hayward, Irwin & Bregman 2005). Malinl was indeed 
discovered in HI in the background of the Virgo cluster (Impey 
& Bothun 1989). How many unknown low surface brightness 
galaxies remain hidden in the local universe, unseen by optical 
observations? ALFALFA should provide the complete census of 
such objects. 

In the last two years ALFALFA has covered most of the Virgo 
cluster (Giovanelli et al. 2006), the nearest rich cluster, provid- 
ing ~ 1000 HI sources in 80 deg 2 , perhaps the highest space 
density of HI sources in the whole ALFALFA survey. The anal- 
ysis of this sample will dramatically enhance our understanding 
of the issues outlined above. The Virgo region is especially use- 
ful because of the existence of a very deep optical catalogue of 
galaxies, the Virgo Cluster Catalogue (VCC) by Binggeli et al. 
(1985), that will allow a comparison of the properties of opti- 
cally and radio selected objects in this region. The overlap be- 
tween the ALFALFA survey and Goldmine, a collection of HI 
pointed observations in the direction of Virgo obtained by many 
observers during the years (Gavazzi et al. 2005), also allows a 
comparison of the relative quality of the two databases. These 
are the issues that the present paper on the Virgo cluster, as seen 
by ALFALFA, wishes to address. A forthcoming paper will an- 
alyze Ha observations obtained recently of several hundred HI 
galaxies discovered by ALFALFA in the Virgo cluster and in its 
immediate surroundings (Gavazzi et al. in preparation), address- 
ing the issue of the environmental dependence of galaxy evolu- 
tion (see a review by Boselli & Gavazzi 2006) by comparing the 
rate of transformation of primordial gas into young stars in the 
Virgo cluster with that in less dense environments. 
The structure of the present paper is as follows. The HI and opti- 
cal selected samples are presented in Section 2. The data analysis 
is presented in Section 3, where the reduction of optical images 
taken from the Sloan Digital Sky Survey (SDSS) is discussed 
in detail. Section 4 gives the comparison between the results of 
ALFALFA and pointed observations available in the literature. 
The results are discussed in Section 5. 

2. The sample 

The present analysis is focused on the ALFALFA blind HI sur- 
vey in the region of the Virgo cluster; more precisely in the area 
of intersection between ALFALFA, as available in Sept 2007, 
and the VCC. The ALFALFA survey provides a sample of HI 
selected objects, and the VCC a list of optically selected galax- 
ies. 

2.1. The Hi selected sample 

The HI selected sample consists of 187 ALFALFA sources 
meeting the following criteria: 

i) their right ascensions lie in the interval 12 /l 08 m 30 I < R.A. < 
12M8 m 20* (ca. 182° < R.A. < 192°) (J2000) corresponding 
to the extent in R.A. of the VCC, and their declinations range 
between 08° < 6 < 16° (J2000), i.e., they lie in the declination 



strip that has been fully mapped by ALFALFA. North of 12 deg 
the catalogued sources are taken from Giovanelli et al. (2007); 
for 08" < 6 < 12" the list of sources is not yet published (Kent 
et al. in preparation); 

ii) they have recessional velocities in the range -1000 < cz < 
3000 km s , bracketing the full depth of the Virgo cluster (see 
Gavazzi et al. 1999); 

iii) they have quality Code = 1, i.e., their statistical significance 
is higher than S/N=6.5. (We do not consider Code = 2 objects, 
i.e., sources with S/N < 6.5 that have optical counterparts with 
optical redshifts matching the HI line redshift.) 

iv) they are not included in the circle of projected radius of 1 .0 
deg centered on M87, where ALFALFA is strongly incomplete 
because of loss of spectral sensitivity due to the bright (220 Jy 
at 1415 MHz) continuum source associated with M87; 

v) they lie outside one resolution element (3.5 arcmin) of any 
other strong source at similar redshift. The six excluded sources 
because of this criterion are: 

122140.1 +143621 = VCC0497 conf by VCC0483 
122140.1 +143621 = VCC1673 conf by VCC1676 
124331.5 +113500 = VCC1972 conf by VCC1978 

Since the mean rms of ALFALFA is 2.1 mJy/beam, the sur- 
vey detection limit (for S/N =6.5) is 0.5 Jy/beam km s _1 for typ- 
ical widths of 40 km s , thus, at the mean distance of the Virgo 
cluster (17 Mpc), using the relation: 

M m = S/N x rms xWx dist 2 x 2.36 x 10 5 (M ) (1) 

these sources have Mm > 10 7 5 M©. The Mm limit, however 
depends on the line width, thus on the system's mass and in- 
clination, as detailed in Section 5.1, resulting in 10 75 M < 
M m limit < 10 8 M Q . 

2.2. The optically selected sample 

From the VCC catalogue (Goldmine version, Gavazzi et al. 
2003), restricted to the region in common with ALFALFA we 
have extracted two subsamples of optically selected galaxies 
that are bona-fide Virgo cluster members, whose distances have 
been assigned following the subcluster membership criteria of 
Gavazzi et al. (1999), updated with new redshift measurements. 
The VCC catalogue contains galaxies as faint as 20 mag, but is 
complete to mj < 18 mag, within a limiting surface brightness 
of 25.3 B mag arcsec~ 2 (Binggeli et al. 1985). From this cata- 
logue we have considered: 

Sample A (deep) consisting of 1 1 12 galaxies with nig < 20 mag; 
Sample B (shallow) consisting of 469 galaxies with oig < 17 
ma£l 278 of which are early-type (dE-E-SO-SOafl and 191 are 
late-type (Sa-Irr-BCD) galaxies. The number of galaxies with 
redshifts is 410/469, including 226/278 early-types and 184/191 
late-types. At the mean distance of the Virgo cluster (17 Mpc, or 
M - m = -31.1 mag) mg < 17 mag corresponds to Mb < -14.1 
mag, thus even Sample B includes dwarf systems. 



1 We have preferred not to extend the shallow sample B to m B = 18 
mag because in the bin 17 < m B < 18 mag only 6% of the VCC galaxies 
have an HI counterpart in ALFALFA, while in the bin 16 < m B < 17 
mag 40% are in ALFALFA. 

2 Notice that the morphological class "dE" might be contaminated by 
some transition objects (dE/dlrr) but we find at most 8 gas-rich ones, as 
shown in Section. 3.4. 



3. The data reduction 
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3.1. HI (ALFALFA) 

For each of the HI selected sources we consider the following 
HI parameters taken from ALFALFA (Giovanelli et al. 2007): 

- the (J2000) coordinates (affected by an error that depends on 
S/N, with a mean of « 24 arcsec); 

- the recessional velocity VP; 

- the observed width of the HI line W. This is the width of the 
profile measured at the 50% level of each of the two peaks, 
corrected for instrumental broadening; 

- the total flux under the HI line: S intP (Jy km s~'); 

For sources optically identified with VCC galaxies, 
we convert SintP into the HI mass using: Mhi = 
2.36 x 10 5 x SintP x (dist) 2 (M ), where the distance is 
determined as in Gavazzi et al. (1999), otherwise using 17 
Mpc, and we correct the observed line width W for projection 
effects, using the ratio of the optical minor- to-major axes: 
W c = W/ sin(acos(b /a)). 

We also estimate the HI deficiency parameter following Haynes 
& Giovanelli (1984) as the logarithmic difference between Mhi 
of a reference sample of isolated galaxies and Mhi actually ob- 
served in individual objects: Defni - LogMni re f-LogMHi obs.- 
LogMni ref has been found to be linearly related to the galaxies 
linear diameter d as: LogMni re f - a + bLog(d), where d (in 
kpc) is determined at the 25'' 1 B-band isophote, and a and b are 
weak functions of the Hubble type. We caution that the Haynes 
& Giovanelli (1984) reference sample of isolated galaxies 
included only relatively large (a > 1 arcmin) UGC objects so 
that the Defyi parameter is poorly calibrated for smaller objects, 
making determinations of the HI deficiency for the smallest 
objects uncertain and likely underestimated (Solanes 1996). 
Furthermore, as discussed in Solanes et al. (2001) galaxies in the 
latest Hubble types (Scd-Im-BCD), for which we have adopted 
a and b parameters consistent with those of Sc, are more subject 
to observational biases than higher surface brightness galaxies. 
The reader should be aware that the determinations of the HI 
deficiency for these objects is highly uncertain. 
For all radio sources we compute the total flux divided by the 
line width: E/// c = SintP/ W in Jy/beam; this parameter is an 
indicator of the "visibility" in the ALFALFA survey (see Figj8]l. 



3.2. Optical & HI (Goldmine) 

For each of the VCC galaxies we take from the Goldmine 
database: 

- the (J2000) optical coordinates (accurate to « 30 arcsec); 

- the morphological type (adopted from Binggeli et al. 1985 and 
from Binggeli et al. 1993); 

- the major (a) and minor (b) axes (arcmin) (from Binggeli et 
al. 1985). These quantities are used to correct the observed line 
widths according to the galaxy inclination; 

- the distance (in Mpc) (computed according to the criteria of 
Gavazzi et al. 1999). The adopted distances for each subgroup 
of Virgo are: 17 ± 0.3 Mpc for cluster A, the North and East 
clouds and the Southern extension, 23 ± 0.5 Mpc for cluster B 
and 32 + 0.9 Mpc for the W and M clouds. 

In order to compare the new HI parameters derived by 
ALFALFA with the plethora of earlier pointed observations, col- 
lected and made as homogeneous as possible by Gavazzi et al. 
(2005)for the Goldmine HI database, we also consider: 

- the measured width of the HI line. The Goldmine measurement 
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Fig. 1. The distribution in celestial coordinates of sources in the 
Virgo cluster. Contours from the X-ray ROSAT all sky survey 
(Bohringer et al. 1994) are given by dashed lines. The solid 
line shows a one degree circle centered on M87, within which 
ALFALFA is complete due to the strong continuum of M87. 
Top: ALFALFA selected sources identified with VCC galaxies 
with n%B < 20 (filled circles). Ten non-identified sources la- 
beled as tidally-disrupted clouds by Kent et al. (2007) (empty 
circles) and 21 sources identified with faint galaxies not listed in 
the VCC (filled triangles) are given. Bottom: Optically selected 
galaxies, consisting of 191 late-type galaxies with ra# < 17 in 
the VCC. 139 that have been detected by ALFALFA are plot- 
ted with filled circles, 15 that have been detected by ALFALFA 
with Code > 2 with filled triangles and another 37 detected 
in Goldmine, but not ALFALFA (as their HI mass is Mhi < 
10 7 ' 5 M Q based on Goldmine) with empty circles. 
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Fig. 2. Comparison of the ALFALFA (blind) and Goldmine 
(pointed) 21 -cm line width measurements. The most deviant ob- 
jects are labeled with their VCC names. The line represents the 
one to one relation. 

differs from that of ALFALFA, as it corresponds to the average 
of the line widths measured at the 20% and 50% levels of the 
peak intensity, uncorrected for instrumental broadening; 
- the hydrogen mass M HI = 2.36 x 10 5 x Sint x (dist) 2 (M ) 
where S int is derived from pointed 2 1 cm observations. It should 
be noted that the 2 1 cm data in Goldmine were obtained from a 
variety of observations taken at Arecibo, some of which lasted 
much longer (typically one hour) than the ALFALFA observa- 
tions (~ 48 sec in the beam), thus reaching a sensitivity about 
one order of magnitude better than ALFALFA. These measure- 
ments are however limited to late-type galaxies. 

3.3. SDSS imaging material and reduction 

For the optical comparisons, we obtained g and i band images 
from the Sloan Digital Sky Survey (DR 5) (Adelman et al. 2007) 
that cover the 187 HI selected sources and the optically selected 
late-type galaxies with oig < 17 (Sample B). Only one galaxy 
(VCC1401) was not available the SDSS. Two or more (up to 4) 
SDSS fields were combined for 13 galaxies (VCC66, 89, 167, 
596, 865, 873, 939, 1588, 1727, 1778, 1932, 2058 and 2070) 
that were present in two or more SDSS images. All frames were 
photometrically calibrated using the wavelength dependent pa- 
rameters contained in the calibration tables associated with the 
individual SDSS images, namely the zero point (ad), the extinc- 
tion coefficient (kk) and the airmass. The effective zero point (in 



Fig. 3. Comparison of the ALFALFA (blind) and Goldmine 
(pointed) 21 -cm line intensity measurements of Virgo galaxies. 
The line represents the one to one relation. 

AB magnitudes) was derived using: ZP e fj = aa + kk * airmass. 
(See the SDSS data release 6 for further details on the flux cali- 
bration.) 

The SDSS image analysis was carried out in the IRAF environ- 
ment and relied on the STSDAS package and on GALPHOT 
(developed for IRAF- STSDAS mainly by W. Freudling, J. 
Salzer, and M.P Haynes and adapted by S.Zibetti and L. Cortese. 
See Gavazzi et al. 2001). 

For each frame the sky background was determined as the mean 
number of counts measured in regions of "empty" sky, and it 
was subtracted from the frame. Sky-subtracted frames were in- 
spected individually and the light from superposed or nearby 
stars and galaxies was masked. 

The 2-dimensional light distribution of each galaxy was fit with 
elliptical isophotes, using the task cphot, a modified version of 
the STSDAS isophote package. Starting from a set of initial pa- 
rameters given manually, the fit maintains as free parameters the 
ellipse center, ellipticity and position angle. The ellipse semi- 
major axis is incremented by a fixed fraction of its value at each 
step of the fitting procedure. The routine halts when the surface 

3 IRAF is the Image Analysis and Reduction Facility made avail- 
able to the astronomical community by the National Optical Astronomy 
Observatories, which are operated by AURA, Inc., under contract with 
the U.S. National Science Foundation. STSDAS is distributed by the 
Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under NASA 
contract NAS 5-26555. 
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brightness found in a given annulus equals the sky rms. The fit 
fails to converge for 46 galaxies with very irregular light distri- 
butions or low surface brightness. In these cases we kept fixed 
one or more of the initial parameters. In 14 instances of exceed- 
ingly low brightness no reliable fit was obtained (for ten HI se- 
lected objects and for four optically selected galaxies, namely: 
VCC 585, 1121, 1582, 1884). 

The resulting radial light profiles were fitted with models of the 
galaxy light distribution. The majority, 72%, were best fit with an 
exponential disk law, while 27% were best fit by a combination 
of a disk and a bulge r 1 ^ 4 law and 1% by a pure de Vaucouleurs 
law (de Vaucouleurs 1948). 

Total magnitudes g, i were then obtained by adding to the flux 
measured within the outermost significant isophote the flux ex- 
trapolated to infinity along either the exponential law that fitted 
the outer parts of most galaxies (pure disks and B+D galaxies), 
or the r 1 / 4 law (dV galaxies). The mean statistical uncertainty 
in the determination of the total magnitude is 0.10 mag. The ef- 
fective radius r e (the radius containing half of the total light, in 
arcsec) and the effective surface brightness fi e (the mean sur- 
face brightness within r e in mag arcsec" 2 ) of each galaxy are 
computed only for the i band. The mean statistical uncertain- 
ties of the determination of r e and [i e are 2 arcsec and 0.15 
mag arcsec" 2 , respectively. The statistical uncertanties refer to 
the errors obtained from fitting models to the observed galaxy 
light distributions; we note that they do not correlate with mag- 
nitude or surface brightness, but rather with the presence of com- 
plex structures (e.g. bars or irregularities in the light profiles). 
Besides the statistical uncertainties, however, the photometric 
parameters of faint galaxies are affected by additional uncertain- 
ties that might depend on local fluctuations of the sky brightness. 
In the following analysis all fits are performed with the maxi- 
mum likelihood method using statistical uncertainties. 
Finally we computed another structural parameter: the concen- 
tration index (C31), defined by de Vaucouleurs (1977) as the 
model-independent ratio between the radii that enclose 75% and 
25% of the total i light. 

We converted the total /-band magnitude to the stellar mass us- 
ing: 

hg(M s , ar /M B ) = -0.152+0.5 18 *(g-i)+bgl(see Appendix A2 
in Bell et al. 2003), where / is the /-band luminosity. The typical 
uncertainty of log(M star /I) is ~ 0.1 dex. 



3.4. The radio-optical identifications 

We cross-correlated the 187 HI selected sources with the opti- 
cally selected sample A (deep) to obtain the optical identifica- 
tions of all HI sources matching a simple positional criterion 
(separation<l arcmin). Fig.Q] top panel, provides the locations 
of these objects. We find 156/187 (83%) identifications with 
VCC galaxies (including 8 sources associated with early type 
galaxies - these coincide with the Code 1 sources found by di 
Serego Alighieri et al. 2007). Ten of the remaining 31 sources 
are not associated with optical galaxies and were found to corre- 
spond to tidally-disrupted systems by Kent et al. (2007) and by 
Haynes, Giovanelli and Kent (2007). Among the other 21 (11% 
of the total sample), 15 are associated with faint galaxies that 
are surprisingly not listed in the VCC in spite of being clearly 
detected and measurable on the SDSS material. Only 6 are not 
visible in the SDSS images; we consider their effective surface 
brightness fainter than 25 mag arcsec" 2 . 

We then cross-correlated the 469 galaxies in the optically se- 
lected VCC sample B (shallow) with the ALFALFA positions, 
separately for the 191 late-type galaxies and for the 278 early - 
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Fig. 4. The Tully-Fisher relation as determined with ALFALFA 
line widths and /-band asymptotic (AB) magnitudes, converted 
to I j and corrected for internal extinction, limited to non- 
deficient galaxies with inclination larger than 35 deg. Filled 
symbols represent spiral galaxies, while empty symbols refer to 
Irr/BCDs. The line gives the slope of the /-band Tully-Fisher 
determination by Masters et al. (2006). 

type galaxies. Fig. [T| bottom panel, provides the locations of 
these objects. Only 4 matches are found with the 278 early-type 
galaxies of Sample B: VCC93, 209, 304 and 355. Among the 
late-types galaxies, 139/191 (73%) match an HI source. Thus 
only 52/191 (27%) late-type VCC galaxies in Sample B do not 
match HI selected sources. Of these, 37 lie below the sensitivity 
threshold of ALFALFA, based on Goldmine; the remaining 15 
are not considered by us in spite of being detected, because they 
have Code>2. 

In any case most of these low Mhi galaxies (32/37 or 46/52 
including the code>2) are relatively bright and HI deficient 
(Def m > 0.6). 



4. Analysis 

4.1. Consistency with previous pointed HI measurements 

The comparison of the line widths determined by ALFALFA 
and by the pointed observations collected in the Goldmine 
database is shown in Fig. [2] There is a tendency ALFALFA HI 
line widths to be up to 50% smaller than those of Goldmine for 
galaxies with W < 100 km s . The difference decreases for 
increasing line widths. This discrepancy reflects the difference 
of the two methods used to measure the line widths in the two 
datasets (i.e. average of 20% and 50% of the peak intensity in 
Goldmine vs. 50% level of each of the two peaks in ALFALFA). 
The line measurements taken homogeneously by ALFALFA 
should provide a reliable determination of the distances to 
individual galaxies using the Tully-Fisher relation (see next 
Section). 

Fig. [3] shows the comparison of the HI mass measurements. 
There is no apparent trend with Mhi- 
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Fig. 5. Hi/optical luminosity correlation plotted separately for 
non-deficient (filled dots) and for objects with Defni > 0.6 
(open symbols). The dashed line gives the (line width depen- 
dent) limit of ALFALFA computed for an inclination of 45 deg. 
The diagonal dotted lines represent the maximum likelihood re- 
gressions fitted to the data, independently for HI deficient and 
normal galaxies. 



4.2. The Tully-Fisher relation 

Using the corrected line width W c as measured by ALFALFA 
and the /-band luminosity derived from the total i^B magnitude, 
and restricting to objects with inclination larger than 35 deg, we 
obtain the /-band Tully-Fisher relation for galaxies in the Virgo 
cluster (see Fig. |4). For this specific purpose it was necessary 
to consider a subsample of non-deficient galaxies, i.e. galaxies 
with the HI deficiency parameter Defni < 0.6. This is because 
in the ram-pressure scenario (e.g. Abadi et al. 1999) HI ablation 
proceeds outside-in, depleting first the gas that is less gravita- 
tionally bound, so that the full width of the measured line pro- 
file underestimates the rotational velocity. After allowing for a 
0.50 mag shift for converting i^B (SDSS) into Ij (Johnson) mag- 
nitudes (according to the conversion table in NED), we apply 
the internal extinction corrections as in Giovanelli et al. (1997). 
For log(W c ) > 2 and excluding the Irr/BCD galaxies, the ob- 
tained relation appears consistent in slope and zero point with 
the template Ij Tully-Fisher determined by Masters et al. (2006). 
This indicates that the distances of the individual Virgo galaxies 
adopted in this work are not badly determined. For log(W c ) < 2 
and for Irr/BCD the dispersion of the relation makes it useless as 
a distance indicator. 



5. Discussion 

5.1. The HI to stellar mass relation 

A basic feature of galaxies is the fact that the HI mass increase 
with increasing optical luminosity is flatter than the direct pro- 
portionality (Roberts & Haynes 1994). Fig. [5] shows the rela- 
tion between these two quantities (given for all optically iden- 




Fig. 6. Gas to star mass ratio separately for non-deficient (filled 
dots) and for objects with Defni > 0.6 (open symbols). The 
dashed line represents the limit of the HI survey computed for 
an inclination of 45 deg. 



tified HI selected sources). Because the present analysis is car- 
ried out with galaxies belonging to a rich cluster, the study of 
the dependence of HI mass on luminosity must be corrected 
for the effects of ram-pressure stripping (that causes HI defi- 
ciency: Giovanelli & Haynes 1985) by excluding HI deficient 
objects (Deffi/ > 0.6). This is why galaxies with "normal" HI 
content are plotted separately from the Hi-deficient objects in 
Fig. [5] (We remind the reader that the definition of HI defi- 
ciency adopted here is a diameter-, not an optical-luminosity- 
based relation). The line fitted to the non-deficient objects (us- 
ing the "maximum likelihood" method) shows that the slope of 
the relation with the optical (/-band) luminosity is < 1, with 
log(M HI /M Q ) = (0.563 ± 0.002) x log(I lum) + (3.68 + 0.02). 
After transforming the optical luminosity into the stellar mass 
(using the relation given in Sect. 3.3) the meaning of the above 
finding is emphasized by plotting in Fig. [6] the ratio of the HI 
to stellar mass versus the stellar mass itself. By focusing only 
on the non-deficient objects it is apparent that the most mas- 
sive galaxies, the giant spirals with stellar masses in excess of 
10 10 Mq, have approximately 5% of their stellar mass in gas, 
intermediate mass galaxies have 10% gaseous mass, but in the 
dwarfs the gas content can exceed the mass of the stars by a fac- 
tor up to 10. Molecular hydrogen is not expected to change this 
trend. Using a luminosity/metallicity dependent CO to Hi con- 
version factor, it has been shown that the molecular gas fraction 
decreases with increasing stellar mass and that the MHijMui 
ratio stays approximately around 15% regardless of the morpho- 
logical type (see Fig. 6 in Boselli et al. 2002). The extrapolation 
of this relation to even lower stellar masses, suggests the ex- 
istence of completely gas-dominated objects, similar to the so 
called "dark galaxies". One of the results of ALFALFA is how- 
ever that, at least in the mass range covered by the survey, there 
are practically no dark galaxies in the Virgo cluster (Kent et al. 
2007; Haynes, Giovanelli & Kent 2007). 
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In both Fig. [5] and [6] we plot the position of Malinl (taken from 
Pickering et al. 1997) to show that we are sensitive to Malinl 
type objects, if they existed in the Virgo cluster. 
As pointed out in Section 2.1, the mass limit of any flux-limited 
HI survey depends on the line width, thus on the system's mass. 
Using the virial theorem, assuming a constant mass/light ratio 
of 4.6 (Gavazzi et al. 1996) and the best fit relation between 
R e and the stellar mass (Fig. [9] and Table 1), we derive a rela- 
tion between W (the line width) and the stellar mass. Using eq. 
1, we obtain the sensitivity of ALFALFA in HI mass units as 
log(M HI /M e ) = 4.25 + 0.37 * logM star - 0.15 (holding for non- 
deficient objects), where the last term represents the correction 
for a galaxy mean inclination of 45°. For less inclined galaxies 
the sensitivity increases, approaching log{Mni IM Q ) = 7.5. The 
relation is plotted in Figs [5] and [6] representing the limiting HI 
mass of sources in ALFALFA as a function of the stellar mass. 

5.2. The surface brightness to stellar mass relation 

One interesting issue that the ALFALFA observations of the 
Virgo cluster allow us to address, as outlined in the Introduction, 
is the question whether optically selection criteria strongly bias 
our knowledge of galaxies. For example, low-surface bright- 
ness objects, particularly giant, low-surface brightness, gas-rich 
galaxies such as Malinl, are strongly undersampled in opti- 
cal catalogues, but should be detected in a radio selected sam- 
ple . To explore this issue we exploit the ALFALFA sample 
in conjunction with the excellent (20 mag limiting magnitude 
and ~ 25.3 mag arcsec -2 limiting surface brightness) VCC cata- 
logue, by comparing the optical structural parameters derived for 
the optically selected and for the HI selected samples. We first 
compare in Fig. [7j the scaling law between the optical (/-band) 
effective surface brightness and the stellar mass (for the subsam- 
ple of the ALFALFA or VCC galaxies that we could measure on 
the SDSS images, i.e., all but 14 objects). The data include the 
Hl-selected sample (filled circles & squares) and the optically- 
selected sample of late-type galaxies with oig < 17 (empty cir- 
cles). The two sets of data appear to be consistent one another. 
The maximum likelihood fit to the optically selected sample is 
plotted together with the lines at +2cr. The horizontal line drawn 
at fi = 23.9 mag arcsec -2 is the i band mean limiting surface 
brightness of SDSS that we have determined as lcr of the sky 
via analysis of hundreds of images. The few measurements that 
lie below this line are meaningless, since they are affected by ~ 1 
mag arcsec -2 uncertainty. There are at most two galaxies (VCC 
307 and VCC 905, highlighted in Fig. [7} that have slightly low 
surface brightness for their mass, and four BCD galaxies (i.e. 
VCC 334, VCC 410, VCC 1313 and VCC 1437) at the oppo- 
site side of the relation, with high surface brightness for their 
mass. There are no others that mimic Malinl. The prototype of 
giant low-surface brightness galaxy, Malinl, plotted for compar- 
ison, deserves some caution. The galaxy is plotted twice (con- 
nected with a line): once at the position given by the discoverers 
(Bothun et al. 1987) with ji/ (V) = 25.5 mag arcsec -2 (that we 
transform into /U o (0 = 25.0 mag arcsec -2 ) at 10 11 LQ, the other 
as measured on the SDSS plates. On this material the low sur- 
face brightness extended disk is not detected, instead the galaxy 
appears as an amorphous enhancement of high surface bright- 
ness (/i e ~ 20.5 mag arcsec -2 ), consistent with the inner disk de- 
tected in the HST measurement of Barth (2007). This raises the 
question whether other Malinl type objects are missing in Fig. 
UJ because their low-surface brightness disks are below the lim- 
iting surface brightness of SDSS. These systems would wrongly 
count as high surface brightness systems because we would de- 




log M star (M ) 

Fig. 7. The fi e vs. M star relation for HI selected late-type galaxies 
(filled circles), optically selected late-type galaxies (empty cir- 
cles) and HI selected early-type galaxies (squares). The diagonal 
lines represent the linear fit to the data with ±2<x. Discrepant ob- 
jects are labeled with their VCC names and Malinl is given for 
comparison. The horizontal line drawn at y. — 23.9 mag arcsec -2 
is the i band mean limiting surface brightness of SDSS. 



tect only their bulges or inner disks, possibly affecting the cor- 
relation between fi e and the stellar mass or even making it spu- 
rious. To make sure that this is not the case we inspected in- 
dividually all galaxies on the SDSS and deeper plate material 
(Goldmine) and we concluded that all high surface brightness 
objects plotted in Fig.|7]are genuine high surface brightness spi- 
rals and BCDs, not bulges associated with missed extended disks 
of much lower surface brightness. The BCDs cannot be mistaken 
for Malinl objects because the associated HI components have 
low mass and widths not exceeding 100 km s -1 , opposite to the 
HI rich and massive Malinl. Even if they were perfectly face- 
on, their much bluer color would discriminate BCDs from pos- 
sible bulges of missed Malinl type objects. We conclude that no 
Malinl -type galaxy is found in the Virgo cluster that we could 
have missed optically or by HI selection, in agreement with 
Minchin et al. (2004), who put severe upper limits on the exis- 
tence of massive low surface brightness objects in the local field. 
Moreover we find a definite trend for the /-band surface bright- 
ness to correlate with the system mass. Extremely low surface 
brightness galaxies exist, but they all lie at the low mass limit of 
HI and optical surveys. At the limiting sensitivity of ALFALFA 
at Virgo (M H1 ~ 10 7 5 ~ 8 M ) there are at most 10/185 HI sources 
that could possibly be identified with galaxies of lower surface 
brightness than the SDSS limit (~ 23.9 mag arcsec -2 ). 



5.3. The distribution of mass and surface brightness 

Fig. [8] contains histograms of the distribution of surface bright- 
ness and luminosity, both HI and optical, as derived for optical 
(red) and HI selected samples (blue). It may be surprising that 
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Fig. 8. Histograms of luminosities (top) and surface brightness (bottom) for the Hi-selected (left) and optically-selected (right) 
samples. The distributions refer to the Hl-selected sample (blue) and to the optically-selected sample limited to m B < 17 (solid red) 
and to m B < 15 (dotted red). Note that the bins at /i e = 25 mag arcsec~ 2 represent upper limits (10 HI selected and 6 optically 
selected objects too faint to be measured in /-band). 



there appears to be very little difference between the red and the 
blue histograms. There is indeed some excess of low surface 
brightness objects in the HI selected sample, represented by the 
10 very faint galaxies that could not be measured on the SDSS 
material and are plotted as upper limits in the 25 mag arcsec~ 2 
bin. (Note that the 10 tidally-disrupted systems were not 
counted among the low-surface brightness objects.) However 
these sources amount to less than 5 % of the total number. 
Similarly in the optically-selected sample there are 4 VCC 
galaxies that could not be measured on the SDSS material. 
The difference is surprisingly small when we compare the 
Hl-selected sample with the mg < 17 optically-selected VCC 
sub-sample. To show that this agreement is not at all obvious we 
compute the distribution also for a shallower subsample of the 
VCC, limited to mg < 15 (red dotted histogram). This time the 
optical (luminosity and surface brightness) distributions differ 
significantly from the Hl-selected distributions, owing to the 
existence of the surface brightness vs. luminosity correlation 
that was investigated in Section 5.2 (see Hig. [7J. Virgo appears 
just at the right distance from us to make a survey with the 
sensitivity of ALFALFA to match almost exactly the VCC 
galaxies of late-type with mg < 17. 



5.4. Other optical structural parameters 

Inspired by the analysis carried out by Cortese et al. (2007) 
who studied the relations between various structural parameters 
in their Fig. 6, we derive additional optical (/-band) structural 
parameters for the HI selected objects and for the nig < 17 opti- 
cally selected late-type VCC sub-sample. These are the galaxy 
"scale" (stellar mass, R e , luminosity or combinations) and 
"form"(jU e , g - i, C31) (following the terminology of Whitmore 
1984). They are plotted in Fig. [9] for sources identified with 
galaxies selected in the blind ALFALFA HI survey (blue) 
and for the optically selected late-type systems from the VCC 
catalogue, limited to mg < 17 (red). The maximum likelihood 
regressions fitted to the data are plotted in Fig.|9]and reported in 
Table 1. 

It is evident from Fig. [9] that the scaling relations well known 
for optically selected galaxies (e.g. color vs. mass, color vs. 
surface brightness, C31 vs. mass; Gavazzi et al. 1996, Scodeggio 
et al. 2002), including the Kormendy relation (R e vs. ju e ), also 
hold true for HI selected objects. We emphasize that this result 
applies to galaxy members of a rich cluster such as Virgo. 
For example the /i e vs. M star relation clearly shows that low 
surface brightness galaxies exist (and are easy to miss), but 
they lie mostly at the low end of the luminosity function, as 
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discussed above. Notice in the color vs. stellar mass and color 
vs. surface brightness panels that the deviant points (faint-red 
objects) are generally dE that have structural and photometric 
properties intermediate between quiescent and star forming 
dwarfs (Boselli et al. 2008) 

Furthermore the frequency distributions of the "scale" and 
"form" parameters are indistinguishable in HI and optically 
selected surveys. For example the frequency of pure disks, that 
have C31 < 3 exceeds the frequency of bulge+disk systems 
(C31 >> 3) in both the HI- and optically-selected samples. 
But galaxies with C31 >> 3 galaxies exist in both samples. 
We conclude that the scaling laws that exist among the optical 
(/-band) structural parameters of disk galaxies do not differ 
significantly in HI and optically selected samples, although 
we should emphasize again that this result is limited to a rich 
nearby cluster. 



6. Conclusions 

By surveying the Virgo cluster, ALFALFA has provided definite 
evidence that: i) HI inhabits galaxies that are structurally similar 
to ordinary late-type galaxies; ii) their HI content can be pre- 
dicted from their optical luminosity; iii) low surface brightness 
galaxies have low optical luminosity and contain small quantities 
of neutral hydrogen; iv) Malinl type galaxies are comfortably 
rare objects (less than 0.5 %); v) there are no "dark-galaxies" 
with HI masses M HI > 10 7 5 ~ 8 M G (see Kent et al. 2007 and 
Haynes, Giovanelli and Kent 2007); vi) less than 1% of early- 
type galaxies contain neutral hydrogen with Mhi > 10 7 5 ~ 8 M 
(di Serego Alighieri et al. 2007). 

Once ALFALFA covers a substantial fraction of its final area, 
including the Virgo cluster in its full extent, a determination of 
the HI mass function will be obtained with unprecedented sig- 
nificance. We expect that the observed HI mass function of the 
cluster can be predicted with sufficient accuracy from the optical 
luminosity function of the late-type galaxies alone, accounting 
for the HI vs. optical luminosity relation discussed in Section 5.1 
(Fig. 5), computed separately for normal and HI deficient objects 
(Kent et al. in preparation). If this is the case we expect that the 
faint-end slope will not differ significantly from that of the opti- 
cal luminosity function of late-type galaxies, i.e. flatter than the 
halo mass function predicted theoretically from the CDM cos- 
mology. 
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Fig. 9. Optical structural parameters for the HI- (blue) and optically-selected mj < 17 late-type galaxies (red). See text for explana- 
tion of parameters. Lines represent the maximum likelihood regressions fit to the data (see Table 1). 



Meyer M. J., et al., 2004, MNRAS, 350, 1195 
Minchin, R.F., et al., 2003, MNRAS, 346, 787 
Minchin, R.F., et al., 2004, MNRAS, 355, 1303 

Pickering, T.E., Impey, CD., van Gorkom, J.H., & Bothun, G.D., 1997, AJ, 
114,1858 

Press, W. H., & Schechter, P., 1974, ApJ., 187, 425 
Roberts, M.S., & Haynes, M.P., 1994, ARAA, 32, 115 
Rosenberg, J. L., Schneider, S.E., 2002, ApJ, 567, 247 

Sabatini, S., Davies, J., Scaramella, R., Smith, R., Baes, M., Linder, S. M., 
Roberts, S., & Testa, V., 2003, MNRAS, 341, 981 

Scodeggio, M., Gavazzi, G, Franzetti, P., Boselli, A., Zibetti, S., Pierini, D., 
2002, A&A, 384, 812 

di Serego Alighieri, S. et al., 2007, A&A (in press) 
Springob, C. M., Giovanelli, R., Haynes, M., 2005, ApJ, 621, 215 
Warren, B.E., Jerjen, H., Koribalski, B.S., 2006, AJ, 131, 2056 
Whitmore, B.C., 1984, ApJ, 278, 61 



Zwaan, M. A., Briggs, F. Ft., & Sprayberry, D., 2001, MNRAS, 327, 1249 
Zwaan, M.A., et. al., 2003, AJ, 125, 2842 
Zwaan, M.A., et. al., 2005, MNRAS, 359, L30 



